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1. INTRODUCTION 

A century-old history of calculation of asymptotics for solutions to Painleve 
equations (usually denoted Pj, j=l,2,...,6) as their variable x tends to infin- 
ity was started by pioneer works by Painleve, Gambler and Boutroux [|I|]. 
In 1980-1981 papers by Jimbo, Miwa and Flashka, Newell initiated the 
so-called Isomonodromy Deformation Method (IDM) which afterwards has 
given remarkable progress in this direction for Pi, P2, and P3. Successful 
attempts to apply the Whitham-like method connected with use of different- 
scaled variables were accomplished in papers 0. This way yields leading term 
of the asymptotic series in form of "modulated" elliptic function (Weierstrass 
p-function for Pi or Jacobi sine for P2) such that parameters determining 
it depend on x in force of certain algebraic or differential equation. At last, 
the Whitham approach in its "classical" interpretation was applied to Pi by 
Novikov, Dubrovin, Moore, Krichever [Q. These authors obtained certain 
ODE usually referred to as "Whitham "or "modulation" equations. In paper 
[^] solutions of Discrete Painleve- 1 eq. were investigated from this point of 
view. The main goal of this paper is to generalize ideas of the latter method 
and expand them to the remaining Painleve eqs. 

2. IDEOLOGY OF THE METHOD 
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Painleve eqs are known to be integrable nonlinear ODE's. Probably, 
the most explicit illustration of this fact is existence of commutative matrix 
representation 



dxLj - d,Aj + [Lj, Aj] =0, j = 1, 2, .., 6, (1) 

where dx = Lj = Lj(u,u' , x, X), Aj = Aj(u,u\ x, \) are 2x2 matrices 
rationally depending on spectral parameter A so that j-th Painleve eq. u" — 
Pj{u,u',x) = is equivalent to (|l]). The matrices Lj, Aj are listed, for 
example, in 0. 

Introduce now variable X and change all variables x explicitly entering 
Lj, Aj to X : Lj = Lj{u,u' , X, X), Aj = Aj{u,u' , X, X). 

Proposition 1. Let e be some real positive number. Then system 

edxLj-dxA, + [L^,Aj]=0 (2) 
is equivalent to system 

dxX = e , s 

u" -Pjiu,u',X) = o 

Proposition 2. Solution of eq. while e = and X = const is 

Uo{x) = fj{T + <i>;'a) (4) 

where r = xU; fj are periodic functions which can be represented in terms 
of Weierstrass or Jacobi elliptic functions (see 0) for any of the six Painleve 
eqs. U = U{X), $ = ^{X);~a{X) are parameters determining the elliptic 
function fj. 

Suppose now e to be small and positive. We shall look for solutions of 
as formal series in e : 

u{x) = Uo{x) + eui{x) + ... (5) 

so that parameters determining the elliptic functions uq = fj obey some 
special ODE called "Whitham" or "modulation" system. Thus, we look for 
the leading term of (|) in form Uo(r, X) = /, + HX); ^{X)) , where 

dxS = U. The Whitham system can be easily derived from eq. @ by simple 
change Udr + €dx and averaging in r : 
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dx det Aj = a22dxlu + audxhi - a^dxhi - a2idxh2 (6) 

where Aj = (amn), Lj = {Imn), m,n = 1,2; bar means averaging over period 
of the elhptic function (p. 

Proposition 3. There exists only unique coefficient in the polynomial detAj 
whose dynamics in X is non-trivial. Denote this coefficient Fj. Therefore the 
Whitham system can be written out as a single first-order ODE on Fj. 

Corollary. The simplest way to find elliptic ansatz is to solve equations 

Fj = consti^ X = const2 (7) 
Proposition 4. The Whitham system induces formal condition 

m'o - Pj{uo,UQ,x) = 0(e)[] 

Obviously, the scaling transformation x ex induces d^X = e — > d^X = 
1 in d^), which provides the following proposition: 

Proposition 5. The function uq specified by (H) and (||) yields leading term 
of asymptotic series for solution of appropriate Painleve equation as x tends 
to infinity. 

3. CONCRETE CALCULATIONS. 

Pi : u" - 3n2 - X = 

det v4i = + 4XA - Fi, where Fi = {u'f - 2u^ - 2Xu; dx det Aj = 
2m + 4A ^ 

dxFi = —2u = —2r]/uj = 2ei + 2{e^ — Ci)^, where E are complete 
elliptic integrals: 

K = S E = J J^^^dz, e = ei, 62 63 are 

roots of the polynomial R^i^z) = Az^ — g2Z — g^, g2 = —X, g^ = —\Fi. 
uo = 2p{x + $; 5(2, 5-3)- 

P2 : u" - 2^3 - ux = Og 

det A2 = 16A^ + 8XA2 + - 4F2, where F2 = (u'f - - Xu^; 

^This condition makes rigorous sense only for smooth function u. However, the method 
describes singular solutions as well. 

^Case with additional constant can be investigated in the same way 
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dx det ^2 = 8A2 + 2X+4m2 ^ dxF2 = -u^ = 2^-e„ = x[\- ^f) , 
where = ^^(X t ^X'' - 4^2) - 1; 

Mo = =F Va;2 - 4^2)5^ (^^^-a; ± Va;^ - 4F2" + $; /c^ , sn is Ja- 

cobi sine. 



P3 : m" + ^m' + sin m = 

det A^ = ^ + + 1, where F3 = ^ - cos m; 

dxF, = -|(F3 + cosm) = ^ (^2 _ 1 + 2f ) , A;2 = (^-Fs ± ^Ff-l)'; 

Mo = -2z log (/csn (if + $; A;)) 

P4 : u" = ^ + ^ + 4:xu'^ + u{a - 4d + 2x^) - f^, where a, c? are 
arbitrary parameters, det = — ^ — XA"' — (d — a + X^)A^ — F4 — p-, 
^4 = ^ - ^ - - M_(X2 + I - 2rf) + X(2d - a) 
9xdet = -A^-2XA2+M2+2XM+a-2c/ ^ 9x^4 = -M2-2XM-a+2d 
and Mo = ±U[C,{Ux + $1) - C(Ux + $2)] ± UC{^2 - ^1) - X, parameters 
U, $2 ~ '^'i can be derived from (|^). The Whitham eq. looks too large to be 
written out explicitly. 



^5 : u" = 2^ -V+f^-r- /"'SttT' is a parameter; 



2MM'(/i - 4) + MV(Ai - 8) 



2m(m-1) V"" / X ' l^x 1^ 2(u-l)' 

det A, = -4X2-^, F5 = [iu' f 

dx det = -8X + where z = ^ 9x^5 = -Az; 

Mo = ^ [p(f/x + $1) + piUx + $2)] 

+ fSfS'u'1) mx + $1) - aux + $2)] + 1^, 

parameters U, $1 — $2, V can be derived from (0). 
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M" 



I u{u—l){u—x) 
a:2(a;-l)2 



1(1+1+^) _ (1 + 1 + ^) ^' 

2 V" u—x J ^ ' \x x—1 u—x J 

y 

(u—x)^ _ 



a, P, 7, 6 are parameters; 
A [(A - 1) (A - x)f det Aq = kik2\^ + FeX^ + SiX + So; 
One can determine and the elliptic ansatz from the following con- 
straint: 
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{u'f x^{x - if + u'ao + (^1 - k2fu^ + 03!*^ + a2U^ + aiu + 9lx^ = 0, 
where: 

\ao = u^x{x-l){ki + k2)+ux [x^iOo + 9i) + x{92 - 9i) - 9o - 92]-9ox^ + 
9ox\ 

- 4^2(^1 - k2){x + 1) + (3A;2 - ki){x{9o + ^1) + ^0 + ^2] - 2^6, 
a2 = a;' [(^0 + ^1) (^0 + ^1 + 4A;i - 4A;2) + 4A;2(2A;2 - + 2x[10A;2 - 6kik2 - 
2kl - 4ki9o 

+^2^1 + 2Fq] + 4fc2(2A;2 - 2ki) + (^0 + ^2)(^o + ^2 + 4fci - 4^2) + 4^6, 

|ai =a;2[(^o+^i)(2A;2-2A;i-^o)+2A;2(A;i-2A;2)]+x[(^o+^2)(2A;2-2A;i-^o) + 

+2A;2(A;i - 2^:2)] - 2^5; 

A;i + A;2 = -(^0 + ^1 + ^2), A;i - A;2 = ^00, « = |(^oo - 1)', P = -i^^, 
7 = 1^?, 5=1(1 + ^1); 

dxFe = -2A;ifc2+Z2(fci-fc2)+^2(^2-fc2)+ - ^^^^f^Sr^^ , t^^^ere 

^0 = ^{M(M-l)(M-a;)2;2+[6'i(M - x) + x92{u - 1) - 2A;2(m - - a;)] z+ 
+k^{u - X - 1) - A;2(^i + x92)}, 
^1 = - (^^{^(^ - - ^)^'+ 

[(6*1 + 6'oo)(it- a;) +2^612(1*- 1) - 2^2(1* - l){u-x)] z+kl{u-x)-k2{9i + 
x92) - kik2}, 

^2 = i^ii^M^ - - ^)^'+ 

-x)+ x{92 + ^oo)(« - 1) - 2A;2(m - l){u - x)] z+kj{u-l)-k2{9i + 
x92) — xkik2}, 

^ u'x{x—l)—8o(u—l){u—x)—6iu{u—x)—02u{u—l) 

2u{u—l){u—x) 
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